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ABSTRACT 

We present a complete solution to the cyclotron-synchrotron radiation 
due to an isotropic distribution of electrons moving in a magnetic field. 
We make no approximations in the calculations other than artificially 
broadening the harmonics by a small amount in order to facilitate the 
numerics. In contrast to previous calculations, we sum over all relevant 
harmonics and integrate over all particle and observer angles relative to 
the magnetic field. We present emission spectra for electron temperatures 
T = 5 x 10 s K, 10 9 K, 2 x 10 9 K to 3.2 x 10 10 K, and provide simple 
fitting formulae which give a fairly accurate representation of the detailed 
results. For T > 3.2 x 10 10 K, the spectrum is represented well by the 
asymptotic synchrotron formula, which is obtained by assuming that the 
radiating electrons have Lorentz factors large compared to unity. We give 
an improved fitting formula also for this asymptotic case. 



1. Introduction 

Radiation by quasi-relativistic and fully relativistic electrons in magnetic fields is 
very common in astrophysics. Examples of sources include supernova remnants (e.g. 
Anderson, Keohane & Rudnick 1995), jets and lobes in radio galaxies (e.g. Carilli et 
al. 1991), hot accretion flows onto neutron stars and black holes (e.g. Narayan, Yi 
& Mahadevan 1995), and relativistic fireballs in gamma-ray bursts (e.g. Paczyhski & 
Rhoads 1993, Meszaros, Laguna & Rees 1993). In all these examples we have thermal 
or non-thermal electrons radiating cyclo-synchrotron radiation in a magnetic field 
which is probably of near-equipartition strength. It is clearly important to be able to 
calculate the radiative luminosities and spectra of these systems. 

Surprisingly, a full solution to this problem does not appear to have been published 
so far, especially for thermal sources. The majority of previous work (e.g. Schott 1912, 
Schwinger 1949, Oster 1960, 1961, Pacholczyk 1970) is limited either to the cyclotron 
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or the synchrotron limit where analytical formulae may be used. These calculations 
are not relevant in the transition zone between cyclotron and synchrotron radiation 
when the velocities of the electrons are quasirelativistic. Petrosian (1981) has obtained 
some approximate formulae in this intermediate regime, but a complete analysis of 
the emission from thermal sources with temperatures in the range 10 9 — 10 10 K can 
only be done numerically. Takahara & Tsuruta (1982) and Melia (1994) have reported 
numerical calculations in the mildly relativistic regime, but both calculations involve 
some approximations and do not agree with each other. Therefore, there is need for 
a more accurate treatment to resolve this difference. In particular, we find significant 
differences between the results quoted by Melia (1994) and the exact calculations 
presented here. 

We present in this paper the complete solution to the cyclo-synchrotron problem 
for a thermal plasma. In we give a qualitative as well as quantitative description 
of the problem, and briefly discuss previous work in this subject. We also present 
an outline of how we solve the problem exactly. Then in §|3] we present the results. 
We calculate the emission due to a thermal distribution of electrons for a range of 
temperatures, and provide analytic fitting functions for the spectra. 



2. Theory 



2.1. Qualitative Features 



A non-relativistic particle moving in a magnetic field radiates with a simple dipole 
pattern (e.g. Rybicki & Lightman 1979) so that an observer sees a sinusoidal electric 
field. The observed spectrum consists of a single delta function at the orbital frequency, 
uj , of the electron. This is the limit of extreme cyclotron radiation. If the particle 
is speeded up, the dipole radiation pattern gets distorted and the observed electric 
field is no longer purely sinusoidal. The spectrum then picks up additional Fourier 
components, or harmonics, which are integral multiples of the fundamental frequency. 
The situation is illustrated in Fig. la. At the same time, u decreases with increasing 
velocity v, according to 

u b eB 
u = — , to b = , (1) 

7 m e c 



where ui b is the cyclotron frequency of the particle, and 7 = 1/ Jl — v 2 / c 2 = 1/ \/\ — (3 2 
is the Lorentz factor. 

With increasing 7, the emission is beamed more and more towards the direction 
of motion of the particle. For large 7, the electric field appears as a series of delta 
function-like pulses, each of width ~ l/^ 2 u bl with succeeding pulses separated by a 
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time 27T7 The Fourier transform of these delta functions gives closely spaced delta 
functions in frequency which approximate a nearly continuous spectrum extending 
up to uj ~ 7 2 w&, and cutting off exponentially at higher u. This corresponds to the 
synchrotron limit, where the spectrum has the universal shape shown in Fig. lb. 

The cyclotron and synchrotron limits have very different spectra, as Fig. 1 
shows. How does the spectrum evolve from one to the other as the particle velocity is 
increased? Answering this question requires detailed calculations and is the topic of 
this paper. 



2.2. Previous Work 



We begin by reviewing a number of well-known results. The power 
(energy/time/steradian/frequency) emitted by an electron moving with a velocity 
parameter between (3 and (3 + d@, in a frequency range du, and at an observer angle 9 
with respect to the magnetic field, is given by (e.g. Schott 1912, Rosner 1958, Bekefi 
1966) 

cos 8 — (3\ 
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(3± sin 9, 



y = moj — uj{1 — (3\\ cos9), 



(2) 

(3) 
(4) 



S(y) is the Dirac delta function, J m (x) is the Bessel function of order m, J' m (x) is its 
derivative, /3m = j3 cos 9 P is the velocity parameter parallel to the magnetic field B, 



where the particle moves at an angle 9 P to the local field, and /3j 
velocity parameter perpendicular to the magnetic field. 



j3 sin 9 P is the 



Each integer m in the summation in equation (fj) corresponds to a harmonic, and 
the presence of the 5-function implies that the emission occurs at discrete frequencies. 
To calculate the power in successive harmonics at low (3 (the cyclotron limit), we 
expand the Bessel functions to lowest order in x. After integrating over all observer 
angles, the emission in each harmonic is given by (Schwinger 1949, Rosner 1956, Bekefi 
1966) 

t 2e 2 u 2 b (m + l)(m^), 



Vn 



-P 1 



(5) 



c (2m +1)! 

This result, which is valid so long as m(3 <C 1, shows that the spectrum decreases 
rapidly with increasing m. Fig. la shows an example with /3 = 0.1, which corresponds 
to this limit. 
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In the opposite limit of 7 ^> 1, which is the synchrotron limit, we are interested in 
large m's, and the Bessel functions can be approximated by modified Bessel functions 
(Bekefi 1966). We then have the familiar result (e.g. Schwinger 1949, Pacholczyk 1970, 
Rybicki & Lightman 1979) 



dE V3e 3 B sin dp 
du 2irm e c 2 



F(X), (6) 



where 



F(X) is given by 



X = — , lo c = - "f 2 uj b sin 6 P , (7) 



F(X)=X J Ks&dt, 



and Ks(£) is the modified Bessel function. By taking appropriate limits of F(X), we 
obtain 

47T /X\ 1/3 

F{x) ^ vfF(l) ( 2 ) ■ x « < 9 > 

F(X) -> ^y^e^X 1 / 2 , X>1. (10) 

Fig. lb shows the synchrotron spectrum for a particle with 7 = 10. 

The above results are for electrons of a given velocity. For a thermal distribution 
we need to integrate the emission over a Maxwellian. This integral can be done in the 
synchrotron limit to give (Pacholczyk 1970), 

e u du = C j— I I -r^- I du ergs s' 1 Hz" 1 , (11) 
K 2 {l/9 e ) \sm9pj 

where 

v 3 7rc 3 6^ m e c^ 

i^fa) is the modified Bessel function, and I(xm) is defined by 

1 °° 

J(x A f) = / z 2 exp(-z) F(x M / z 2 )dz. (13) 

In writing equation fllTD , we have included the proper normalization for a relativistic 
Maxwellian distribution; in the limit of high temperatures, this normalization is 
equivalent to that used by Pacholczyk (1970). The limiting behavior of I(xm) for small 
and large Xm are straightforward. The former has been worked out by Pacholczyk 
(see Eq. fl£^)), and the latter by Petrosian (1981) (see Eq. ([X4]). Equation ([□]) 
corresponds to particles at a fixed angle 9 P to the magnetic field. For an isotropic 
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particle distribution, we need to perform another integral over d(cos 9 P ). The asymptotic 
dependencies of the result for small and large % are again easy to calculate and are 
discussed in Appendix A. 

For the mildly relativistic case, Petrosian (1981) has obtained the following analytic 
approximations for the emission spectrum from a thermal distribution of electrons as a 
function of the observer angle 9: 



e w (6) du^^C T ^ *,„ , exp (-6.75 1/3 a^ 3 ) du, x &e > 1, 



and 



6Va K 2 (l/6 e 



e u {6) du - &I*2*I 2 t?C 
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e e K 2 (i/e e 



(xd< 



v3/2 



(14) 



1 + cos 2 6»" 
sin 2 



exp 



-Xln 



ex# e sin 9 



du, 9 e < 1. 



(15) 



These results are valid only with the additional condition u/ub 3> 1, the regime 
Petrosian (1981) was interested in. For xm 3> 1, equation ( |i~4"| ) is exactly the same as 
equation ([TT| ) (cf. Eq. |A4j) , and we see that the approximations made by Petrosian 
(1981) in this regime are the same as the extreme synchrotron approximations. 

The extreme synchrotron results described above provide quite an accurate 
representation of the exact results for electron temperatures ^ 3 x 10 10 K, while the 
series result given in equation (5) is very good at temperatures ^ 10 8 K. For 10 8 K 
^ T ^ 3 x 10 10 K, we cannot use either of these limiting results but must solve the 
complete cyclo-synchrotron problem. 



2.3. The Complete Solution 



We start by writing the integrals that need to be carried out. To determine the 
luminosity, = dE/du, due to a particle moving with velocity parameter /3, we must 
integrate T]u(fl,9) (Eq. (|)) over observer angles, and a given particle distribution: 



1 2?r 1 2tt 1 

dE 2 



= — = — J d(3n((3) J dcPpJ d(cos9 p ) J d<p J d(cos9) Vw (0, 9), (16) 

-1 

where ^($,6) is given by Eq. (0). The subscript p refers to the particle, and n(/3) is 
the velocity distribution of the particles, which is taken to be isotropic. The factor of 2 
in front is because we integrate over only half the range of cos9 p , and the factor of l/Air 
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comes from the angular normalization of the (isotropic) particle distribution function. 
The integrals over <j) p and <j) are trivial, giving 

111 
' 2tt f d(3n{f3) f d{cos9 p ) f d(cos9)r] LJ 0,9). (17) 



du 



'V 

oo-i 



For a fixed velocity parameter j3 and frequency u, we numerically evaluate the two 
innermost integrals in Eq. (|lj]), as well as the sum over harmonics in the expression (2) 
for 77aj(/3, 9). We repeat this calculation for various values of f3 and u and tabulate the 
results. The results can then be convolved with any isotropic velocity distribution to 
obtain the spectrum L w . In this paper, we restrict ourselves to a relativistic Maxwellian 
n((3) and present detailed results for this particular case. 

The evaluation of Eq. flTT| ) involves a delta function (Eq. (|2|)), which determines 
the precise frequencies at which radiation is observed. For a given harmonic m, the 
delta function implies that there is emission only at 

muj 

" = (l-/3||eosfl) - (18) 

Since m takes on only integer values, this means that for a given {3 and cos#, emission is 
observed only at a discrete set of uj. Unfortunately, the discrete nature of the emission 
poses a serious numerical difficulty since it requires infinite resolution in frequency. In 
order to make the numerics tractable, we replace the delta function with a smooth 
broadening function, f(uj), which is nonzero over a finite frequency range u c ± Auj. 
Here uj c is the central frequency where we wish to evaluate the emission, and Auj is a 
broadening width, which we adjust. With the smoothing function, the delta function 
in equation (2) becomes 

S(y) = 6 [mu - uj{1 - (3 {l cos9)] = - ± f( x ), (19) 

L J Ub 1 — P\\ cos 9 

where 

X=- = — ■ (20) 
For /(x), we choose the functional form 

15 



fix) 



WA X 



A X 2 ^ ™ \A X 



(21) 



which has the property that f(ui c ± Au) = f'(ou c ± Auj) = 0. 



In our calculations we set A\ = a X , and we have found that a choice a = 0.05 
allows the delta function to be broadened sufficiently to stabilize the numerics without 
losing too many details of the harmonic structure of the emission. For each value of Xc 
of interest, we calculate the power emitted into a width 2A X centered at Xc , counting 
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only those harmonics which fall within this width. With Eqs. (0), and (19), Eq. (|T 
thus becomes 

'*W\ } Ja ro, } Jt ns},, n, fix') 



2vr (1^1 J d{3n{{3) J d{cos6 p ) J rf(cos( 



X 





oo / nnB 

12 



1 — (3\\ cos 6* 

ft, cin ft \ 

(22) 



m _-i v sin^ y \ 1 — /3y cos ^ y \ 1 — /?|| cos 6 1 

For a Maxwellian velocity distribution with temperature T , we have 

9 

w(7) rf 7 = kTK ^g j /^7 2 exp(-7/0 e )rf7, (23) 

or 

2 

nC9) = fcT ^ 2 e (1/ge) /?V exp(- 7 / W, (24) 

with 

/OO /*1 
n ( 7 ) rf 7 = J n ((3) d(3 = 1. (25) 

Here /c is Boltzmann's constant, # e = kT/m e c 2 , and i^Oc) is the modified Bessel 
function. 

The procedure to calculate the full cyclo-synchrotron emission spectrum is as 
follows. Using Eq. (0), we fix the values of (3, cos9 p , cos 9, and sum over the 
harmonics responsible for emission over the frequency range x ~ Ay to x + A\> This 
gives the total emission into the given observer angle 9 from particles at angle 9 p . We 
then integrate over all observer (cos 9) directions, to get the emission into all of observer 
space, and integrate over all particle (cos# p ) directions to obtain the total emission per 
particle at the given value of x due to an isotropic distribution of particles moving with 
the given j3 or Lorentz 7. We repeat the calculation for various values of 7 and x to 
obtain a table of values. 

The present calculations differ from previous work by Tsuruta & Takahara (1982) 
and Melia (1994) in that we integrate over all observer angles whereas those authors 
restricted themselves to a particular observer angle. 



3. Results 



3.1. Transition from Cyclotron to Synchrotron. 



Figure 2 shows plots of the emission as a function of scaled frequency x f° r various 
particle velocities, after the angular integrations in Eq. fl22|) have been performed. For 
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comparison, the dashed lines show the synchrotron formula as given in Eq. (|6|). We 
see that the exact numerical results deviate considerably from the limiting synchrotron 
formula in panels (a)-(c), where the harmonic cyclotron-like emission is quite evident. 
The harmonics are broadened, partly because of our broadening function f(x) (Eq. 
(|2lf) , and partly due to Doppler broadening induced by (3\\ (see Eq. (fl8|)). The former 
effect dominates in panel (a), but the latter is more important in panel (c). 

By panel (d), which corresponds to (3 = 0.9, 7 = 2.3, the harmonics have merged 
completely at higher values of \ and the synchrotron approximation is becoming quite 
good. Yet higher values of 7 (panels (e)-(f)) make this trend more apparent. However, 
at low values of \ the harmonics continue to be present and there are deviations from 
the synchrotron formula. 

Figs. 3a-d show polar plots of the emission as a function of the observer angle 9. 
In these figures, the particle has 7 = 10.3 and is moving in a helix at an angle 9 P to the 
magnetic field with cos^p = 0.3. The field points along the positive y-axis. Fig. 3d 
shows the emission pattern that is observed at a large frequency x = 100- According to 
Fig. 2f, which corresponds to this value of 7, the emission at this \ is practically equal 
to the synchrotron formula. The reason is clear from Fig. 3d. We see that the emission 
is tightly beamed along the particle direction, and the beam has a half angle ~ I/7, 
exactly as assumed when deriving the synchrotron formula. Varying cos# p does not 
change this behavior and, therefore, on integrating over cos 9 P we find a total emission 
which agrees very well with the synchrotron formula. 

In strong contrast is the case shown in Fig. 3a, which corresponds to the same 
values of 7 = 10.3 and cos6 p = 0.3, but has a much lower value of x = 0-8- The 
emission pattern now is not strongly beamed about the particle's velocity. Instead 
we can clearly see a set of harmonics, each with a width determined by f{x)- As 
cos 9p varies, the individual beams get longer or shorter, but the harmonic character 
remains. Therefore, the integrated emission over all cos6 p deviates significantly from 
the synchrotron approximation (as Fig. 2f shows at log(x) — —0.1). Figs. 3b, c, show 
how the transition in the beaming occurs as x increases. 



3.2. Emission by a Thermal Distribution of Electrons 

The primary output from our calculations is a table of r] x values as a function of 
electron velocity (3 (or equivalently 7) and dimensionless frequency x- This table can 
be convolved with any isotropic velocity distribution to calculate the corresponding 
spectrum. In the following we describe the results for a thermal Maxwellian distribution 



- 9- 



3.2.1. The Ultra-relativistic Regime 



We first consider the ultra-relativistic regime, where kT 3> m e c 2 . Pacholczyk 
(1970) derived the formulas given by Eqs. ([11]) and QTB]), for the emission from a 
relativistic Maxwellian distribution of electrons with a fixed particle angle 8 P . For an 
isotropic distribution of particles, we define a new function I'(xm) 



and substitute I'(xm) for I(xm) in Eq. (JTTJ) . 

Pacholczyk (1970) has numerically calculated I(xm) over a wide range of xm and 
tabulated the values. In addition he has shown that in the limit of small xm 

'^-^M^T- 2 - bWxr - (27) 

In the opposite limit of % ^> 1, Petrosian (1981) finds 

I(x M ) -> 2.5651 exp(-l. 8899 xff). (28) 
We can combine these two limiting extremes into the following simple fitting function, 



1 92 9977\ 1 /c 

I(x M ) = 2.5651 ( 1 + — ^ + -— - exp(-1.8899a# 3 ), (29) 

M X M / 



where we have optimized the coefficient 1.92 in the middle term so as to minimize 
the error. This fitting function has a maximum error of 0.39%, which occurs when 
xm ~ 63. Figures 4a, b, compare the fitting function with the exact numerical values 
of I{xm) and show the residuals. 



In the case of I'(xm), we show in Appendix |X] that as xm — > 

I'(x M ) -> 2.1532 x M 2/3 , (30) 

and as oo 

I'(x M ) -> 4.0505 x M 1/& exp(-1.8899xlf ). (31) 
Once again we combine these two limits to obtain a fitting function: 

Tl . , 4.0505 / 0.40 0.5316\ / 1/3 x . , 

I'(x M ) = —tjtr 1 + + —jj^r exp (-1.8899 4f) • (32) 

This function has a maximum error of 2.7% at xm ~ 160. Figs. 4c, d, compare the 
fitting function with the numerical values and show the residuals. 
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3.2.2. The Mildly Relativistic Regime 



Using our tabulated values of Tjy^ CIS Bb function of (3 and Xi we nave computed 
emission spectra for isotropic thermal distributions of electrons with temperatures in 
the range, 5 x 10 8 K < T < 3.2 x 10 10 K. We see that at the highest temperatures the 
emission is very similar to I'(xm), equation (|32"D , but there are significant deviations at 
lower temperatures, especially at small x- We have obtained a set of fitting functions 
corresponding to each of the temperatures for which we have calculated the spectrum. 
Each function is of the form 

4.0505a / 0.40/? 0.5316 7 ^ / 1 ____ i/ 3 \ 
M{x M ) = — 1 + —[jr + 172— exp (-1.8896 x^ J , (33) 

X M V X M X M ' 

where a, /3, 7, are all adjustable parameters which we have optimized so as to minimize 
the square of the deviation of M(xm) from the numerically calculated results. Table p] 
shows the optimized parameters we obtained at the various temperatures. We expect 
that as T — > 00, a, j3, and 7 should all — > 1, since M(xm) must approach I'(xm)- 
Indeed we see that this is the case in Table [I]. At lower values of T, however, the 
parameters are very different from 1. This is to a large extent because the fitting 
function is attempting to fit the harmonic "bumps" in the spectrum 

In terms of the fitting function M(xm), the optically thin cyclo-synchrotron 
emission from a thermal plasma at temperature T is given by (cf. equation fllTD), 

e^du = C R(x, T) du ergs s" 1 Hz" 1 , (34) 

with 

RlV 77 = 

K 2 (l/6 e 



xM{x M ) , _v 

R(X,T) = r- , , (35) 



and the power is given by 

w£u> = XL X = uj b CxR(x,T) ergs s _1 . (36) 
For temperatures T > 3 x 10 10 K, M(xm) should be replaced by I'(xm) given in 



equation (|32|), which is equivalent to setting the fitting constants a, ft, 7 equal to unity. 



Figs. 5, 6, show the fits at various temperatures along with the residuals, and 
Table shows where the maximum errors occur, both over the range log(x) < 1 and for 
log(x) > 1- We see that the errors are particularly severe for log(x) < 1 because of the 
harmonic oscillations which are impossible to fit in detail with a simple function such 
as (|33|) . However, in many applications, the synchrotron emission will be self-absorbed, 
and one would be interested primarily in log(x) > 1- We see that the errors here are 
much less severe. 



We have compared our detailed numerical results with those of Petrosian (1981), 
Takahara & Tsuruta (1982) and Melia (1994). We find good agreement with the former 
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two papers. The small differences in our results can be explained by the fact that they 
only considered a single observer direction whereas we have averaged over all directions. 
We do, however, find a serious discrepancy with Melia's results. For instance, our 
calculation shows that for T = 10 10 K, the spectrum peaks at log(x) ~ 2, whereas the 
calculation by Melia indicates that the peak is at log(x) ~ 1 and that the emission 
is insignificant at log(x) — 2. In fact, Melia's results seem to indicate that the cutoff 
frequency for the synchrotron emission is essentially independent of the temperature, 
whereas it is clear from basic principles that the cutoff must increase rapidly with 
increasing temperature. 



4. Conclusion 

The detailed calculations described in this paper bridge the gap between the 
limits of nonrelativistic cyclotron emission and ultrarelativistic synchrotron emission. 
We have calculated the emission due to an isotropic distribution of charged particles 
moving in a magnetic field, and have shown how the spectrum changes as a function 
of the particle Lorentz factor 7 and the dimensionless frequency \ (defined in Eqs. 20 
and 1). Figure 2 shows some results for selected cases. Included in our calculations are 
all the details of the harmonic emission. This is important at low values of 7, and for 
low frequencies even at high values of 7. Also, we isotropically average over observer 
directions relative to the magnetic field. 

Having calculated emission spectra for an array of values of 7, we have integrated 
the spectra over an isotropic relativistic Maxwellian distribution of particle velocities 
to calculate the spectrum of cyclo-synchrotron emission due to a thermal plasma in 
a magnetic field. The results are shown in Figs. 5 and 6 for temperatures ranging 
from 5 x 10 8 K to 3.2 x 10 10 K. We have obtained fitting functions M[xm) with three 
fitting constants, a, /3, 7 (see Eqs. 33-35 and Table 1), which provide a fairly accurate 
representation of the numerical results. These fitting functions allow the spectrum to 
be calculated with reasonable accuracy for any temperature > 5 x 10 8 K. The errors 
decrease with increasing temperature in the manner indicated in Table 2. 

The thermal cyclo-synchrotron spectra presented in this paper agree with most 
previous results, except that our calculations are more complete since the spectra 
have been integrated over all particle and observer angles. At highly relativistic 
temperatures, our results agree with those given by Pacholczyk (1970), while at mildly 
relativistic temperatures and for frequencies u ^> u>b, °ur results agree with those of 
Petrosian (1981) and Takahara & Tsuruta (1982). There are minor deviations in the 
results which can be traced to the fact that the previous authors specified a fixed 
direction of the observer relative to the field rather than averaging over all directions. 
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Our results do, however, differ significantly from the calculations presented by Melia 
(1994) and we have been unable to understand the reason for the discrepancy. 

Finally, we note that the basic output of our calculations is a table of emission 
spectra for isotropic particles of fixed velocity (3 or Lorentz factor 7. In this paper we 
have concentrate on one application of this table, namely the calculation of thermal 
cyclo-synchrotron spectra from thermal plasmas with Maxwellian velocity distributions. 
The tabulated results could be convolved with any other isotropic electron distribution 
function, e.g. a power law distribution, to calculate the corresponding spectrum. 
Our work thus provides a "ready to use" table for determining the cyclo-synchrotron 
emission from any astrophysical source with an isotropic particle and magnetic field 
distribution. 
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Barbara). RN thanks the ITP for hospitality. 



A. Asymptotic Formulae for I(x) and I'(x) 



A.l. I{x) 
The definition of I(x) is given in Eq. (fT3|): 

oo 

I(x) = - / z 2 exp(-z) F(x/^ 2 )^. (Al) 
o 

By using the approximation of F(x) given in Eq. (|9|), Pacholczyk (1970) showed that, 

J(x)^--^Jf-) - 2.5593a;- 2 / 3 , x < 1. (A2) 

x 9^3 V2/ 



In the opposite limit of x 3> 1, we use the approximation of F(x) given in Eq. (JlOf) to 
write 

oo 

Jfr) ^ - ^ 1/2 / exp[-(^ + x/z 2 )] dz. (A3) 
V 2 x 7 2 
o 

Employing the method of steepest descent, Petrosian (1981) showed that this integral 
can be evaluated to give 

2tt 



I{x) - ^exp[-(2 1 / 3 + 2- 2 / 3 )x 1 / 3 ], 

-> 2.5651 exp(-1.8899x 1/3 ), x > 1. (A4) 



A. 2. /'(x) 

The function J'(x) is defined by 



We first consider the case for x <C 1, where we have J(x) — > 2.56x~ 2//3 . Setting 
x — > x/ sin^, we obtain 

J ' (3;) ^ 2 -|/S Sin2/3 ^ (cOS ^' (A6) 
o 

which can be evaluated to give 

~ 2.153x" 2/3 , x<l. 
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For x >> 1, we have I(x) — > 2.5651 exp(— 1.8899a; 1 / 3 ), and setting x — > x/ sin6>, we 
obtain 

i 

/'(a) = 2.5651 y exp(-1.8899x 1/3 /sin 1/3 ^)rf(cos^). (A8) 
o 

We now use the fact that most of the emission comes from 9 ~ ir/2. Setting 9 = + 7r/2 
we obtain 

i 

I\x) -> 2.5651 y exp(-1.8899x 1/3 /cos 1/3 0)rf(sin0). (A9) 
o 

Setting cos0 = (1 — y 2 ) 1 ^ 2 , and Taylor-expanding up to y 2 , we finally obtain 

i 

I'(x) -> 2.5651 y exp(-1.8899a; 1/3 ) exp(-1.8899a; 1/3 y 2 /6)^, (A10) 



2.5651 exp(- 1.8899a; 1 / 3 ) / — vr — / 6 



2 V 1.8899 V a; 1/3 ' 

4.05047^ exp(- 1.8899a; 1 / 3 ), x > 1. 



Table 1: Optimal values of the parameters for different temperatures. 



1 (A J 




1 [X) 






a 


O 

P 


7 


5 x 10 


0.0431 


10.44 


16.61 


1 X ru 


1 1 91 


-1U.00 


Q 1 P.Q 


2 x 10 9 


1.180 


-4.008 


1.559 


4 x 10 9 


1.045 


-0.1897 


0.0595 


8 x 10 9 


0.9774 


1.160 


0.2641 


1.6 x 10 10 


0.9768 


1.095 


0.8332 


3.2 x 10 10 


0.9788 


1.021 


1.031 



Table 2: List of errors and where they occur for log(x) < 1 and for log(x) > 1, for 
different temperatures. 



T (K) 


1 / \ ^ i 

log(x) < 1 


% Error 


1 / \ ^ i 

log(x) > 1 


% Error 




iO&(.Xmax .error ) 




'■O&yXmax. error ) 




k w 1 n8 

e> x 1U 


U.l 


/i /in 
44U 


O 

z.o 


or -i 
OO.l 


1 x 10 9 


0.05 


33.6 


3.3 


16.6 


2 x 10 9 


0.0 


5.9 


3.0 


5.3 


4 x 10 9 


0.33 


6.7 


2.8 


2.4 


8 x 10 9 


0.33 


5.3 


1.9 


0.9 


1.6 x 10 10 


0.0 


7.0 


1.4 


0.56 


3.2 x 10 10 


0.0 


7.4 


1.6 


0.55 
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Figure Captions 

Figure 1: (a) Cyclotron emission for a particle with j3 = 0.1. The emission is normalized 
to the total emission in all harmonics, (b) Synchrotron emission for a particle with 
7 = 10 from standard synchrotron theory. 

Figure 2: Plots of scaled emission log(7 x ) against log(x) for (a) (3 = 0.2, (b) (3 = 0.3, 
(c) (3 = 0.6, (d)P = 0.9, (e) 7 = 5.3, (f) 7 = 10.3, (g) 7 = 40.3, (h) 7 = 100.3. The 
vertical axis is the scaled emission, l x after the angular integrations in Eq. ( p2|) have 
been performed. l x = L x (c/e 2 u%) = (c/e 2 ^) (cf. Eq.(|22"D). The dashed lines show 
the synchrotron limit given by Eqn. (||). 

Figure 3: Polar plot of emission from a particle moving with a Lorentz factor 7 = 10 
at an angle of cos 6^ = 0.3. The magnetic field is oriented along the positive y-axis and 
cos# = corresponds to the x-axis. The values of x are ( a ) X — 0.8, (b) x — 1-4, (c) 
X = 5.0, (d) x = 100.0. 

Figure 4: (a) Plot of the analytic and numerical values for I(x) and (b) the 
corresponding residuals, (c) Plot of the analytic and numerical values for I'(x) and (d) 
the corresponding residuals. 

Figure 5: (a) Plots of Eqn. ( |3"4| ) (dashed lines) and the numerical calculation (solid 
lines) for (a) T = 3.2 x 10 10 , (b) T = 1.6 x 10 10 , (e) T = 8 x 10 9 , (f ) T = 4 X 10 10 , and 
their corresponding percentage errors (c), (d), (g), and (h). 

Figure 6: (a) Plots of Eqn. ( |3~4"D (dotted lines) and the numerical calculation (solid 
lines) for T = 2 x 10 9 , (b) T = 1 x 10 9 , (e) T = 5 x 10 8 , and their corresponding 
percentage errors (c), (d), and (f). 
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